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The human substantia nigra (SN) is thought to consist of two functionally distinct 
neuronal populations — dopaminergic (DA) neurons in the pars compacta subregion and 
GABA-ergic neurons in the pars reticulata subregion. However, a functional dissociation 
between these neuronal populations has not previously been demonstrated in the awake 
human. Here we obtained microelectrode recordings from the SN of patients undergoing 
deep brain stimulation (DBS) surgery for Parkinson's disease as they performed a 
two-alternative reinforcement learning task. Following positive feedback presentation, we 
found that putative DA and GABA neurons demonstrated distinct temporal dynamics. DA 
neurons demonstrated phasic increases in activity (250-500 ms post-feedback) whereas 
putative GABA neurons demonstrated more delayed and sustained increases in activity 
(500-1000 ms post-feedback). These results provide the first electrophysiological evidence 
for a functional dissociation between DA and GABA neurons in the human SN. We discuss 
possible functions for these neuronal responses based on previous findings in human and 
animal studies. 

Keywords: substantia nigra, human, dopamine, GABA, neuron, reinforcement learning 



1. INTRODUCTION 

Animal studies have shown that the substantia nigra (SN) 
consists of two functionally distinct neuronal populations — 
dopaminergic (DA) neurons in the pars compacta subregion and 
GABA-ergic neurons in the pars reticulata subregion. DA neurons 
have been shown to encode reward prediction errors with pha- 
sic bursts of firing, that occur when there is a mismatch between 
obtained and expected outcomes (Schultz et al, 1997; Bayer and 
Glimcher, 2005). These DA bursts are thought to guide rein- 
forcement learning by adjusting synaptic strength in downstream 
regions following unexpected outcomes (Reynolds et al, 2001; 
Tsai et al., 2009). In contrast, GABA neurons are involved in 
inhibitory regulation of various brain structures including frontal 
cortical regions (via the thalamus), premotor brainstem nuclei 
and midbrain DA neurons (Carpenter et al., 1976; Hikosaka and 
Wurtz, 1983;TepperetaL, 1995; Hennyetal, 2012). Despite these 
advances in the animal, the functional role of human SN neurons 
has not been elucidated. 

Patients undergoing deep brain stimulation (DBS) surgery 
for the treatment of Parkinson's Disease offer a rare opportu- 
nity to directly study the functional properties of human SN 
neurons (Jaggi et al., 2004). Two previous studies in patients 



undergoing DBS suggest a functional role for the human SN in 
reinforcement learning. First, it has been shown that a subset of 
neurons in the SN demonstrate phasic bursts of activity follow- 
ing unexpected rewards, consistent with a reward prediction error 
(Zaghloul et al., 2009). Second, microstimulation applied in the 
SN following rewards alters learning by enhancing the reinforce- 
ment of preceding actions (Ramayya et al, 2014). In both studies, 
the observed learning-related neural and behavioral patterns were 
presumed to reflect the function of a healthy subpopulation of 
DA neurons in the region. Although histochemical studies have 
shown that DA and GABA neurons co-exist in the human SN 
(Damier et al., 1999a), a functional dissociation between these SN 
neural populations has not previously been demonstrated. 

In this study, we sought to directly compare the response pro- 
files of DA and GABA neurons recorded from the human SN so 
as to assess whether these neuron groups represent functionally 
distinct subpopulations. We obtained recordings from 25 sub- 
jects as they performed a two-alternative reinforcement learning 
task where they selected between stimuli that carried distinct 
reward probabilities and received positive or negative feedback 
following each choice. We extracted neuronal spiking activity 
from each unit and identified putative DA and GABA neurons 
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based on the physiological properties of their recorded waveforms 
(Joshua et al, 2009; Matsumoto and Hikosaka, 2009; Ungless 
and Grace, 2012). If DA and GABA neurons demonstrate dis- 
tinct task-related responses, it would suggest that they represent 
functionally distinct neuronal populations. 

2. MATERIALS AND METHODS 

2.1. ELECTROPHYSIOLOGICAL RECORDINGS 

We obtained intra-operative microelectrode recordings from 25 
Parkinsonian patients undergoing surgery for the implantation of 
a deep brain stimulator (DBS) in the subthalamic nucleus (STN). 
Patients who volunteered to take part in the study provided their 
informed consent during preoperative consultation and received 
no financial compensation for their participation. Per routine 
clinical protocol, Parkinson's medications were stopped on the 
night before surgery ( 12 h preoperatively); hence subjects engaged 
in the study while in an OFF state. The study was conducted in 
accordance with a University of Pennsylvania Institutional Review 
Board-approved protocol. During surgery, intra-operative micro- 
electrode recordings (obtained from a 1 |xm diameter tungsten tip 
electrode advanced with a power-assisted microdrive) were used 
to identify the substantia nigra (SN) and the STN as per routine 
clinical protocol. We obtained microelectrode recordings sampled 
at 25 kHz using a StimPilot recording system (16 bit analog-to- 
digital converter) and Spike2 data acquisition software (Figure 1; 
targeting and recording details are reported elsewhere; Moyer 
et al., 2007). In this study, we present data captured from the SN 
as subjects performed the reinforcement learning task described 
below (see "Reinforcement learning task"). 

2.2. REINFORCEMENT LEARNING TASK 

Subjects performed a two-alternative probability learning task 
which has been previously used to study reinforcement learning 
and value-based decision making (Figure 2; Frank et al, 2004, 
2007; Zaghloul et al., 2012). During the task, three pairs of sym- 
bols (denoted here by letters: AB, CD, EF) were presented in 



random order, and subjects were instructed to choose one of the 
two stimuli on each trial (Figure 2B). Selections were made by 
pressing buttons on handheld controllers placed in each hand. 
The three stimulus pairs were characterized by different relative 
rates of reward (AB, 80% vs. 20%; CD, 70% vs. 30%; EF, 60% 
vs. 40%). Reward rates associated with each symbol were deter- 
mined randomly prior to each session and were fixed throughout 
the experiment. Probabilistic feedback followed each choice. In 
the event of positive feedback, the screen turned green, and the 
sound of a cash register was presented. In the event of negative 
feedback, the selection screen turned red, and an error tone was 
presented. Each trial consisted of presentation of the stimuli, sub- 
jects response, and a 2 s display of feedback. Subjects were asked to 
make selections which maximized their probability of obtaining 
positive feedback. As in previous reinforcement learning studies 
in the human SN (Zaghloul et al, 2009; Ramayya et al, 2014), 
there was no monetary payout and the provided feedback was 
virtual. 

The rationale for including three item pairs with distinct rel- 
ative reward rates is two-fold. First, we wanted to encourage 
learning throughout the session. Second, it allowed for the study 
of subthalamic nucleus neurons during decision conflict in a 
subsequent experiment. When possible, subjects first performed 
the task during the preoperative consultation, but in all cases, 
the task was reviewed with subjects on the morning of surgery. 
Further instructions were provided prior to beginning the task 
intra-operatively. During surgery, subjects performed the task on 
a laptop placed comfortably in front of them while the microelec- 
trode was positioned in the SN. We aligned behavioral data with 
neural recordings by sending sync-pulses to the neural record- 
ing system from the behavioral laptop as participants performed 
the task. Some participants were bilaterally implanted with DBS 
electrodes and performed two intra-operative sessions of the task. 
The 25 subjects performed 32 sessions in total with a mean (±SD) 
of 123 (±7.1) trials per session. Each session typically lasted ~ 
15 minutes based on participants' response times. 





Waveform duration (ms) 



FIGURE 1 | Intra-operative electrophysiological methods. (A) During deep 
brain stimulation (DBS) surgery, a microelectrode is advanced into the 
substantia nigra (SN) to identify the inferior border of the subthalamic nucleus 
(STN). Top-an example pre-operative MRI scan (sagittal view) overlaid with a 
standard brain atlas and estimated microelectrode position is shown. This 
figure is adapted from Zaghloul et ai. (2009). Bottom-an example 500 ms 
band-pass (400-3000 Hz) filtered voltage trace is shown. We extracted 



neuronal spiking activity from each microelectrode recording by identifying 
spikes in the filtered signal that demonstrated sufficient separation from 
background noise (Materials and Methods). (B) We identified putative DA (n 
= 13, dark gray) and GABA (n = 10, light gray) units based on their baseline 
firing rate and waveform durations. Left: mean waveforms from DA and 
GABA fast-spiking units. Width represents standard error of mean (s.e.m). 
Units that did not fall in either category are marked with open circles. 
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2.3. EXTRACTING NEURONAL SPIKING FROM MICROELECTRODE 
RECORDINGS 

From each microelectrode recording, we extracted neuronal 
activity using the WaveClus software package (Quiroga et al., 
2005). We band-pass filtered each voltage recording from 400 to 
5000 Hz and manually removed periods of motion artifact. We 
identified spike events as positive or negative deflections in the 
voltage trace that crossed a threshold that was manually defined 
for each recording (~4 SD about the mean amplitude of the fil- 
tered signal). The minimum duration between consecutive spike 
events (censor period) was set to be 1.5 ms. Spike events were sub- 
sequently clustered into units based on the first three principal 
components of the waveform. Noise clusters from motion arti- 
fact or power line contamination were manually invalidated. To 
ensure neuronal isolation, we filtered units based on established 
measures of isolation quality (Isol; Neymotin et al., 2011). We 
rejected units if greater than 0.025 of their inter-spike intervals 
were refractory period violations (<3 ms) or if units were poorly 
separable from background noise in feature space (Isolde < 4). If 
multiple units on a channel met the aforementioned criteria, but 
were poorly separated from each other (IsoImv < 4), they were 
considered together as a multi-unit, which is appropriate for our 
analyses because DA and GABA neurons are typically regionally 
clustered in the SN (Henny et al., 2012). We identified a total of 
42 units. Seven units were excluded because of poor separation 
from background noise and/or refractory period violations. Of 
the remaining 35 units, two units were poorly separated from each 
other and were combined into a multi-unit. Thus, our dataset 
consisted of 33 single-units (Isolde = 6.31 ± 1.13; mean ± SD), 
and one multi-unit (IsoIgG = 7.53, mean IsoIjvn = 2.78). These 
data were identified from 17 of the 25 subjects; 18 sessions yielded 
one unit, whereas 8 sessions yielded two units. 

2.4. IDENTIFYING PUTATIVE DA AND GABA ACTIVITY 

To understand the function of SN DA and GABA neurons, we 
sought to extract the activity of these neuronal populations from 
microelectrode recordings. Because pars compacta and pars retic- 
ulata are largely interspersed in the primate SN (Poirier et al., 
1983), the location of the microelectrode relative to any anatom- 
ical landmarks is typically not used to isolate activity from these 
neuronal populations (also, see Menke et al, 2010). Instead, non- 
human primate electrophysiology studies usually identify puta- 
tive DA and GABA units based on the properties of extracellular 
spike waveforms recorded on the microelectrode (Fiorillo et al., 
2013). Previous studies which have combined electrophysiologi- 
cal recordings with pharmacological manipulations (Schultz and 
Romo, 1987) or histochemical techniques (Henny et al, 2012) 
have shown that DA neurons exhibit slow firing rates and broad 
waveforms, whereas GABA neurons display fast firing rates and 
narrow waveforms (Ungless and Grace, 2012). From each unit, 
we estimated baseline firing rate by computing the mean firing 
rate over the entire recording session and waveform duration by 
measuring the peak-to-trough duration (Barto et al., 2004). We 
identified putative DA units as those which displayed baseline fir- 
ing rates slower than 15 Hz and waveform durations >0.8 ms, and 
GABA units as those which displayed baseline firing rates faster 
than 15 Hz and waveform durations <0.8 ms; similar parameters 



have been used in a prior non-human primate study (Matsumoto 
and Hikosaka, 2009). For the multi-unit in our dataset, we con- 
sidered the baseline firing rate to be the average baseline firing rate 
of the two contributing units to account for the artificial elevation 
in firing rate that results from combining units. 

For each DA and GABA unit, we computed smoothed fir- 
ing rates during each trial by convolving the spike train with 
a Gaussian kernel (half-width = 75 ms). To aggregate firing 
rate responses across units, we computed normalized firing rate 
responses for each unit. Specifically, we computed a distribution 
of mean firing rates shown by the unit across all trials (0-1000 ms 
post-stimulus and — 500-2000 ms surrounding response). We 
z-scored the smoothed firing rate during each trial based on 
the mean and standard deviation of this distribution. The time 
intervals used for the normalization process rarely overlapped 
because subjects demonstrated a mean reaction time of 2047 ms 
(±855 ms). 

2.5. STATISTICAL METHODS 

For all statistical analyses, we aggregated activity within each unit 
and studied changes in firing rate across units. We studied fir- 
ing rates from each unit in non-overlapping 250 ms windows 
(0-750 ms following stimulus presentation, and — 500-1500 ms 
surrounding response trials), that were chosen a priori based 
on prior animal (Schultz et al., 1997; Cohen et al, 2012; 
Pan et al, 2013) and human (Zaghloul et al, 2009) stud- 
ies of midbrain DA and GABA activity. To assess whether 
DA and GABA units demonstrated distinct temporal dynam- 
ics, we performed a 2 x 2 ANOVA following the three task 
events (stimulus presentation, responses resulting in positive 
and negative feedback). We considered time interval and neu- 
ron type to be fixed effects. To account for variability that 
may result from obtaining multiple samples from each pop- 
ulation, we included neuron number as random effect nested 
within the neuron-type fixed effect. We performed post-hoc f-tests 
to identify specific changes in neural activity, and corrected 
for multiple comparisons using a false-discovery rate (FDR) 
procedure. 

3. RESULTS 

We obtained microelectrode recordings from the substantia nigra 
(SN) of 25 patients (16 males, mean age = 57.36) undergoing 
deep brain stimulation surgery for the treatment of Parkinson's 
disease (PD). As per routine clinical procedure, microelectrodes 
were advanced into the substantia nigra (SN) in order to identify 
the inferior border of the subthalamic nucleus, the target for the 
stimulating electrode (Figure 1A; Jaggi et al., 2004; Zaghloul et al., 
2009). From each SN recording, we extracted neuronal spiking 
activity and identified putative DA (n = 13, mean rate = 4.56 Hz, 
mean duration = 0.87 ms) and GABA (n = 10, mean rate = 
25.0 Hz, mean duration = 0.62 ms) units based on their base- 
line firing rates and waveform durations (Materials and Methods, 
Figure IB). 

As we obtained recordings, subjects performed a two- 
alternative probability learning task where they were asked to 
select between pairs of Japanese characters by pressing but- 
tons on hand-held controllers. Immediately following each 
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response, they probabilistically received positive or negative feed- 
back (Figure 2A). Each stimulus carried a distinct probability 
of reward and each pair always consisted of a high-probability 
and a low-probability stimulus. During each session, subjects 
were presented with three item pairs that varied in their relative 
reward rates (80/20, 70/30, and 60/40). Subjects were instructed 
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80/20 70/30 60/40 

t Relative reward rate 

FIGURE 2 | (A) Reinforcement learning tasks. During surgery, subjects 
performed a two-alternative reinforcement learning task where they were 
asked to select between pairs of Japanese characters by pressing buttons 
on hand-held controllers. Immediately following each response, positive 
feedback (green screen, sound of cash register) or negative feedback (red 
screen, error tone) was probabilistically provided. An example positive and 
negative feedback trial are illustrated. Subjects were informed that each 
stimulus carried a distinct probability of reward and that their goal was to 
maximize positive feedback during the session. (B) Behavioral performance. 
During each session, subjects were presented with three stimulus pairs 
that varied in their relative reward rates (80/20, 70/30, and 60/40). To index 
subjects' learning during the task, we measured their bias toward selecting 
the high probability item during the final 10 trials of a given pair. Subjects 
reliably demonstrated a bias on the 80/20 pair (0.69) and a modest bias on 
the 70/30 pair (0.6). We did not observe a bias on the 60/40 pair (0.55). 
Error bars represent s.e.m across subjects. See main text for statistics. 
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FIGURE 3 | Distinct responses from DA and GABA units following 
positive feedback. We studied aggregate normalized firing rates from DA 
in = 13, dark gray) and GABA (n = 10, light gray) units in relation to three task 
events — stimulus presentation (A) subject responses that resulted in positive 
feedback (B), and negative feedback (C). We observed distinct responses 



to select stimuli that maximized their probability of receiving 
positive feedback. To index learning on a particular item pair, 
we measured the tendency that subjects demonstrated toward 
selecting the high-probability item during the last 10 presenta- 
tions of that item pair (Figure 2B). We found that subjects reliably 
demonstrated such a tendency during the 80/20 pair [0.69, tyo) = 
4.64,p < 0.001]. Subjects showed a trend toward such a tendency 
on the 70/30 pair (0.60, p = 0.08), but not on the 60/40 item pair 
(0.55, p > 0.2). 

To compare the functional properties of DA and GABA units, 
we studied aggregate normalized firing rates from each popula- 
tion aligned to three task-related events — stimulus presentation, 
responses associated with positive feedback, and responses asso- 
ciated with negative feedback (Figure 3). We separately exam- 
ined neural activity following responses associated with positive 
and negative feedback because DA units have been shown to 
demonstrate opposing responses during these trials (Zaghloul 
et al., 2009). To compare responses from the two groups dur- 
ing these three conditions, we binned firing rates from each 
unit in non-overlapping 250 ms windows (0-750 ms following 
stimulus presentation, and — 500-1500 ms surrounding response 
trials) and applied two-factor ANOVAs with neuron-type and 
time-interval as fixed effects. We included neuron number as 
random effect nested within the neuron-type fixed effect to 
account for variability that occurs when obtaining multiple sam- 
ples from each population (see Statistical Methods). Following 
positive feedback presentation, we observed a significant inter- 
action between neuron-type and time-interval [F(7, 183) = 6.02, 
Mean Squared Error (MSE) = 0.81, p < 0.001] suggesting that 
DA and GABA neurons demonstrated distinct temporal dynamics 
during these trials. Post-hoc f-tests revealed that DA units demon- 
strated greater firing rates than GABA units during the 250- 
500 ms time interval [f( 2 i) = 2.37, p = 0.028] whereas GABA 
units demonstrated greater firing rates than DA units dur- 
ing the 500-750 and 750-1000 ms time intervals [t( 2 i)'s > 



from DA and GABA units following responses associated with positive 
feedback, but not following stimulus presentation or responses associated 
with negative feedback. Firing rate responses were smoothed using a 
Gaussian-kernel (half-width = 75 ms). Width of each response represents 
s.e.m across units. 
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2.52, p's < 0.029; false-discovery rate (FDR) corrected p's < 
0.07]. We did not observe significant interactions between 
neuron-type and time-interval during stimulus presentation or 
following negative feedback (p's > 0.16). Thus, we observed dis- 
tinct responses from DA and GABA neurons following positive 
feedback presentation, but not following stimulus or negative 
feedback presentation. 

To assess whether differences between DA and GABA firing 
rates following positive feedback were driven by changes in DA 
activity, GABA activity or both, we studied changes in each pop- 
ulation's firing rates from baseline. We selected the following 
time intervals of interest based on the results of the previ- 
ous analysis: 250-500 ms ("early," when DA activity was greater 
than GABA activity) and 500-1000 ms ("late," when GABA activ- 
ity was greater than DA activity). For DA units, we observed 
increased firing rate from baseline during the early time interval 
[f(i2j = 2.15, p = 0.052], but did not observe significant changes 
in firing during the late time interval (p > 0.2). For GABA units, 
we observed the opposite pattern — we did not observe signifi- 
cant changes in firing during the early time interval (p > 0.2), 
but observed significant increases in firing rate during the late 
time interval [f(9) = 3.29, p = 0.009]. Thus, the major changes in 
neural activity following positive feedback presentation included 
an early increase in DA activity and a late increase in GABA 
activity. Example DA and GABA units are shown in Figures 4, 5, 
respectively. 

4. DISCUSSION 

We studied neuronal activity in the SN of patients undergoing 
DBS surgery for the treatment of Parkinson's disease as they 
performed a two-alternative reinforcement learning task. During 
each trial of the task, subjects were presented with a pair of stim- 
uli, selected one of the stimuli by pressing buttons on hand-held 
controllers ("response"), and immediately received positive or 
negative audio-visual feedback. We identified putative DA and 
GABA neurons based on the physiological properties of their 
extracellular waveforms, and compared the functional properties 
of the two populations during the task. 

4.1. DA AND GABA NEURONS IN THE HUMAN SN ARE FUNCTIONALLY 
DISTINCT 

Our main finding was that DA and GABA neurons demonstrated 
distinct temporal dynamics following responses that resulted in 
positive feedback. Whereas DA neurons demonstrated phasic 
bursts in activity (250-500 ms post-feedback), GABA neurons 
demonstrated more delayed and sustained increases in activity 
(500-1000 ms post- feedback). These results provide the first elec- 
trophysiological evidence for a functional dissociation between 
DA and GABA neurons in the human SN. Whereas prior histo- 
chemical studies have shown that DA and GABA neurons co-exist 
in the human SN (Damier et al, 1999a), the only direct evidence 
for a functional dissociation between these neural populations has 
come from animal electrophysiology studies (DeLong et al, 1983; 
Schultz et al., 1997). Our findings provide a bridge between these 
studies by demonstrating a functional dissociation between these 
neural populations in the human SN. As such, our results provide 
electrophysiological support for neuro-computational theories of 



human basal ganglia function that ascribe distinct roles to these 
neural populations during learning and decision-making (Bogacz 
and Gurney, 2007). 

4.2. FUNCTIONAL SIGNIFICANCE OF PHASIC DA BURSTS 

Animal electrophysiology studies have shown that DA neurons 
demonstrate phasic bursts of activity that correlate with reward 
prediction errors (Schultz et al., 1997; Bayer and Glimcher, 2005). 
Enhancement of these DA bursts via electrical microstimulation 
(Reynolds et al., 2001) or optogenetics (Tsai et al, 2009) results 
in enhanced learning, suggesting a causal relation between pha- 
sic DA bursts and learning. However, several factors limit the 
generalizability of these studies to human behavior. First, ani- 
mal learning is typically studied following primary rewards and 
punishments (e.g., juice and airpuffs), whereas human learning 
is often motivated by higher-order abstract rewards (e.g., rational 
and social goals). Second, animals in these studies have typically 
undergone long periods of intense training, whereas much of 
human learning occurs in novel situations. 

Recent studies in patients undergoing DBS surgery for 
Parkinson's disease suggest a functional role for phasic DA bursts 
in human reinforcement learning. Zaghloul et al. (2009) demon- 
strated reward prediction error-like responses in a subset of 
SN neurons that electrophysiologically resemble DA neurons 
described in animal studies (putative DA neurons). The current 
study functionally validates the use of these electrophysiological 
criteria by showing that putative DA neurons demonstrate dis- 
tinct post-reward responses from other neurons in the region. 
Consistent with our findings, Ramayya et al. (2014) found that 
microstimulation applied near SN neuronal populations that 
showed post-reward bursts of activity and broad waveforms 
resulted in altered learning. Generally, our finding that puta- 
tive DA neurons demonstrated post-reward bursts in activity 
(Figure 4) is consistent with their hypothesized role in providing 
reinforcement following rewards (Glimcher, 2011). 

Similar to the Zaghloul et al. (2009) study, we observed 
DA bursts 250-500 ms following feedback, which is later than 
DA bursts typically observed in animal studies ( 100-250 ms; 
Niv and Montague, 2008). The more delayed latency might be 
attributed to the presentation of abstract audio-visual rewards, 
rather than primary rewards, each of which might engage DA 
neurons through distinct processes (prefrontal vs. brainstem 
mechanisms, respectively, Glimcher, 2011). Unlike the Zaghloul 
et al. (2009) study, however, we did not observe clear evi- 
dence that post-reward DA bursts represented a reward predic- 
tion error (although, see Supplementary Material). This may be 
because subjects demonstrated limited learning during the task. 
Additionally, whereas Zaghloul et al. (2009) observed DA pauses 
during the 150-300 ms post-feedback interval, we did not observe 
reliable decreases in activity across DA neurons (although, see 
Figure 4B). This discrepancy may be explained by the fact that 
the negative feedback condition in the Zaghloul et al. (2009) study 
was associated with an absence of reward, whereas in our study, it 
was associated with the presentation of a salient negative stimu- 
lus. Previous animal studies have shown that pauses in DA activity 
are less frequently observed following the presentation of aversive, 
salient stimuli (Matsumoto and Hikosaka, 2009). 
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FIGURE 4 | Example DA units. Three representative DA units are shown. 
For each unit, average waveform (top left), inter-spike intervals on a 
logarithmic time scale (bottom left, vertical line indicates 3 ms), smoothed 
rate (half-width = 75 ms) and raster following responses associated with 



positive (middle) and negative feedback (right), respectively, (vertical line 
indicates response). Width of smooth rate represents s.e.m. Baseline firing 
rates, waveform durations for the three units are as follows. (A) 1.31 Hz, 
0.82 ms (B) 6.91 Hz, 0.85 ms (C) 3.35 Hz, 0.92 ms. 



4.3. FUNCTIONAL SIGNIFICANCE OF GABA ACTIVITY 

In contrast to DA neurons, GABA neurons demonstrated delayed, 
and sustained increases in activity following positive feedback. 
These patterns are consistent with findings from animal studies 
that have shown sustained changes in midbrain GABA activity 



following visual stimulus and reward presentation (Handel and 
Glimcher, 2000; Sato and Hikosaka, 2002; Joshua et al, 2009; 
Cohen et al., 2012). We speculate that these post-feedback GABA 
responses are related to a reciprocal interaction with DA neu- 
rons. Previous work has shown that GABA neurons demonstrate 
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FIGURE 5 | Example GABA units. Three representative GABA units are shown. Same conventions as in Figure 4. Baseline firing rates and waveform 
durations are as follows. (A) 25.6 Hz, 0.67 ms (B) 27.0 Hz, 0.75 ms (C) 28.3 Hz, 0.39 ms. 



increased firing rates when exposed to dopamine (Waszczak and 
Walters, 1983), suggesting that DA neurons may exert excitatory 
control of GABA firing. Conversely, SN GABA neurons exhibit 
inhibitory projections onto midbrain DA neurons, and may exert 
inhibitory control over DA neurons (Tepper et al, 1995; Lobb 
et al, 2011; Henny et al, 2012; Pan et al., 2013). Then, follow- 
ing a phasic DA burst, GABA neurons might display an increase 
in firing rate that might act to regulate DA firing and suppress 



subsequent DA phasic bursts. GABA responses might be more 
prominent following positive compared to negative feedback if 
the majority of DA neurons that provide inputs to GABA neurons 
demonstrate preferential increases in phasic activity following 
positive feedback compared to negative feedback. Although the 
majority of SN GABA neurons reside in the pars reticulata, a sub- 
set of GABA neurons are also known to exist in the pars compacta 
region (Nair-Roberts et al., 2008; Ungless and Grace, 2012). 
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Some GABA neurons also demonstrated robust pauses in 
activity soon after feedback was presented (see Figure 5). Pauses 
in GABA-ergic activity typically suggest a release of inhibition 
on downstream structures, and have been classically observed 
during movement and saccade generation (DeLong et al., 1983; 
Hikosaka and Wurtz, 1983). These pauses in activity are thought 
to decrease inhibition on ("disinhibit") downstream motor struc- 
tures (e.g., superior colliculus; Carpenter et al., 1976), and allow 
for the execution of a movement. Thus, the observed GABA 
pauses may be related to some movement expressed by subjects 
immediately following the presentation of salient sensory stim- 
uli during the feedback condition (possibly orienting saccades; 
Hikosaka and Wurtz, 1983). However, we are unable to test this 
hypothesis because we did not monitor eye movements during the 
study. Alternatively, the observed pauses in GABA activity may be 
related to decreased inhibition on DA neurons that would facili- 
tate post-feedback DA bursting (Luscher and Ungless, 2006; Lobb 
etal, 2011). 

4.4. LIMITATIONS 

We note several limitations to our study. First, we are 
unable to provide direct histochemical evidence that these 
electrophysiologically-identified neural subgroups reflect distinct 
neuronal populations. However, there is a large body of evidence 
from animal studies suggesting that these electrophysiological 
criteria may be used to identify distinct midbrain neuronal popu- 
lations (Ungless and Grace, 2012). As such, several animal studies 
rely on electrophysiological criteria alone to identify functional 
subpopulations within the midbrain (Matsumoto and Hikosaka, 
2009; Fiorillo et al., 2013). Second, the population we stud- 
ied in this experiment-patients undergoing DBS for Parkinson's 
disease-is known to have degeneration of neurons in SN. Ideally, 
one would like to study the function of SN neurons in healthy 
human subjects, but at present such recordings may not be 
ethically obtained in any other human population. Converging 
evidence from histochemical (Damier et al., 1999b) and electro- 
physiological studies (Zaghloul et al., 2009; Ramayya et al., 2014) 
in patients with Parkinson's disease and in animals (Hollerman 
and Grace, 1990; Zigmond et al, 1990; Wang et al, 2010) indi- 
cate that a significant population of viable DA neurons remain in 
the Parkinsonian SN. We suggest that the observed DA and GABA 
responses reflect activity from the subpopulation of healthy neu- 
rons that remain in the SN. 
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